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^ ■ ABSTRACT 

["T I ' We solve the problem of propagation and dissipation of Alfvenic turbulence in 

■ a model solar atmosphere consisting of a static photosphere and chromosphere, 

transition region, and open corona and solar wind, using a phenomenological 
>■ ■ model for the turbulent dissipation based on wave reflection. We show that most 

of the dissipation for a given wave-frequency spectrum occurs in the lower corona. 
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CN ' and the overall rms amplitude of the fluctuations evolves in a way consistent 

Q ■ with observations. The frequency spectrum, for a Kolmogorov-like slope, is not 

^ ■ found to change dramatically from the photosphere to the solar wind, however 
it does preserve signatures of transmission throughout the lower atmospheric 

P^' layers, namely oscillations in the spectrum at high frequencies reminiscent of 

O ' the resonances found in the linear case. These may disappear once more realistic 

' couplings for the non-linear terms are introduced, or if time- dependent variability 

ci ■ of the lower atmospheric layer is introduced. 

> : 

^ . Subject headings: MHD - waves - turbulence - Sun: solar wind 

1. INTRODUCTION 

In situ measurement of magnetic and velocity field fluctuations from Helios and Ulysses 
have revealed a broad developed spectrum for frequencies ranging from 10~^ Hz to 10~^ Hz. 
Typically, a strong correlation between m agnetic field and velocity fluctuations in this dis- 



tance range persists (IMangeney et al.lll99ll ) corresponding to an outwardly propagating spec- 



trum. It is well known that nonlinear terms couple Alfven waves propagating in opposite 
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directions. Also, the basic nonlinearity in homogeneous MHD in the presence of a majority 
of one type of waves forces the evolution with time to increase the dominance, preferentially 
dissipating the minority component in a process called dynamical alignment (Veltri et al. 
1980) which is not observed in the solar wind. Therefore the presence of a well-developed 
spectrum together with a preferred direction of propagation has remained a mystery. The 
question that naturally arises therefore concerns the drivers for the continuing and anomalous 
(compared to homogeneous MHD predictions) nonlinear cascade in this outwardly dominant 
case. Among the possible drivers of a nonlinear cascade in the solar atmosphere are com- 
pressible effects, which couple Alfven waves with slow and fast modes, or couplings due to the 
stro ng gradients in the atmosphere. Among the fi rst are phenomena suc h as parametric de- 



cay 



Pruneti fc Velli 1997: Del Zanna et al. 



200ll ) and wave-steepening (ISuzuki fc Inutsuka 



20051 ). Gradients transverse to the mean magnetic field directions lead to phase- mixing, 
i.e. development of small scales in directions perpendicular to that of propagation. Finally, 
the gradients due to stratification cause wave-refiection, which naturally produces the waves 
propagating i n opposite d irecti on required for the classical incompressible cascade, as first 
suggested by I Velli et al.l (119891 ). Disentangling the role of all of these processes at once 
would require fully 3D calculations in a realistic atmosphere model, a feat beyond present 
numerical capabilities. We therefore f ocus here on the role of wave-refiection, which has been 



extensively studied in the linear case ( Heinemann &: OlbertI 



1980; 



Lerov 



among the first) while less so in the nonlinear one (IMatthaeus et al 



1983 



19801: lHollwedll978 



Some constraints on the frequency spectrum and the energies for the outward and inward 
propagating components are derived from the observations. The Alfvenic fiuctuation power 
spectrum in the fast solar wind evolves with distance (R) not self-similarly with a power-law 
dependence on u with slope -1 and -5/3 at low and high frequencies respectively. The two in- 
tervals are separated by a critical frequency (a;*) which depends on R. Identifying the fiuctu- 
ations with Alfven waves it is useful to adopt the Elsasser variables z"^ = u=i=sign(i?o)b/ y/Anp 
(corresponding respectively to outward and inward propagating Alfven waves if the mean 
magnetic field Bq is pointing outward from the sun) The energy per unit mass residing in 
the outward and inward propagating modes {E^ = |z^p res pectively) both decreas e with 
distance and for R < 2.5 AU oc R'^-^^ and R- oc E~°-^^ Jsavassano et al.l boOObI ) . The 
normalized cross helicity, dc = {E^ — E~)/{E~^ + E^) which accounts for the imbalance be- 
tween the outward and inward component, also evolves with distance and it is approximately 
equal to one in t he inner solar wind, it d ecreases for R > 0.4 AU and oscillates around ~ 0.4 
for R > 2.5 AU JSavassano et al.lboOOah . 



It must be recalled that the observed features can not be explained either by linear 
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propagation theory (including reflection) or by MHD turbulence separately. 
A linear analysis applied to the solar wind case shows that low frequency waves {uj < 
experience the strongest reflection in the photosphere, chromosphere and corona (iHoUweg 



Q-^ Hz\ 



19781 . Il98ll : ISimilon &: Zarghanj|l992l ). Their flux at the transition region is greatly reduced 



(even if a considerable power is transmitted to the corona) and in the outer (super sonic) solar 



wind the radial dependence of cTc is similar to the observed at higher frequencies (jVelli et al. 



199l|). 



On the other hand the high frequency waves (10~^ B.z < u < 10^^ Hz) are almost completely 
transmitted (even if in the photosphere and chromosphere their reflection is relatively high), 
both the and E~ energies decrease faster than the scaling observed and finally ~ 1 in 
the outer solar wind (IVelli et al.lll99ll ). 

The dynamics of a well developed turbulent state in the expanding solar wind has been 
studied as well and ordering of the characteristic time scale which should effectively favor 
the development of a turbulent cascade in planes perp endicular to the direction of wave 



propagation (along the magnetic f ield) has been found (IZank et al 



1996; 



Matthaeus et al. 



1998 . 11999I : iDmitruk et all l2001al . I2OO2I : IPmitruk fc Matthaeud l2003l : lOughton et al 



2001 



2OO4J ). Numerical mode ls capable of reproducing the o bserved Uc profiles in the s upersonic 



part of t he solar wind (|z 



1984; Tu 1988; Velli et al 



ou fc Matthaeud Il989l . Il990l ) or the spectral evolution (jTu et al. 



19891 ) necessarily use ad hoc assumption and simplification, and 



even if considerable advances have been made, a complete understanding of the s olar tur- 
bulent spectrum and the solar wind acceleration (ILi et al.lll999l : iHabbal et al.lll995l ) has not 
been achieved. 



Here, we investigate the combined effect of wave reflection and turbulent dissipation in 
order to understand the relative importance of linear and nonlinear effects on the overall 
evolution of the fluctuation amplitudes. Comparison of the numerical results with some ob- 
servations give some constrains on the fields at the photospheric and coronal level for which 
data are still missing with implications for numerical models of solar wind acceleration. 
In the context of a reflection driven turbulent cascade process another interesting issue con- 
cernes the evolution of the turbulent spectrum. If one supposes that the Alfven waves are 
injected at the photospheric base at a well defined frequency or with a given correlation time 
one would expect to find a signature of this characteristic time-scale in the observed spec- 
trum at 1 AU (or in other words, one can ask if discrete modes and turbulence can coexist 
( iDmitruk et al.l l2004j )). No injection frequency is observed in the solar wind spectrum so 
one can ask if both the turbulent evolution and the frequency-dependent transmission prop- 
erties of the solar atmosphere and wind can efficiently smooth this supposed strong forcing 
signature. 
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We integrate numerically the equations for the velocity and magnetic field fiuctuations 
(written in terms of the Elsasser fields) for a stationary model atmosphere with a supras- 
pherically expanding wind, form the photosphere to one AU for a set of frequency chosen 
in the range 10~^ Hz < < 10~^ Hz. Each wave is identified via its frequency while a 
phenomenological nonlinear term is added to the equations in order to account for both 
turbulent dissipation and frequency coupling. 



2. THE MODEL 

The equations describing the propagation of Alfven waves in an inhomogeneous sta- 
tionary medium can be derived from the Magnetohydrodynamic equations (MHD) un- 
der the hypotheses of incompressible adiabatic transverse fluctuations. The velocity (u) 
and magnetic field fiuctuations (b) can be combined to form the Elsasser variables = 
u =p s^gn{'Bo)h/^y4'I^p which describe Alfven waves propagating outward (z"*") or inward 
(z^). Bo stands for the average magnetic field while p is the mass density. In terms of these 
variables the equations for the two fields read: 

^ + [(U ± Va) • V]z± + (z^ ■ V)(U T V.) + 

± ^{z^ - z±)[V ■ V. T ^(V ■ U)] = -(z^ ■ V)z±, (1) 

where U is the mean wind speed and the Alfven speed is Va = Bo/^/47rp, co-linearity be- 
tween magnetic and gravitational field is assumed. On the right hand side we have grouped 
the nonlinear terms (except the total pressure, which in the limit of incompressible fluctua- 
tions can also be written as the product of z+, z and their gradients). In the linear part of 
the eq. [2] we can recognize a propagation term (II) and two terms accounting for reflection 
due to the variation of the properties of the medium, one isotropic (IV) while the other (HI) 
involves variations along the fluctuations' polarization. 




Fig. 1. — From top left, clockwise: wind speed (solid line) and Alfven speed (dotted line), 
numerical density, temperature and expansion factor as function of heliocentric distance for 
the modeled atmosphere. 
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The chromosphere and the photosphere are modeled as a static layer, 2400 km thick, 
with the magnetic field organized in fiux tube in supra-spherical geometry with constant 
temperature. The density varies almost exponentially and the magnetic field varies accord- 
ing to the fiux tube expansion (A) in order to reproduce the properties of a coronal hole in 
the quiet Sun (jHoUweg et al.lll982l ). Across the transition region the density falls off by two 
orders of magnitude, the wind passes from a speed of km s~^ to 8 km s~^ while the mag- 
netic field strength is continuous (about 10 G). The corona also expands supra-spherically 
and its temperature profile is chosen to fit observations (see fig. [T]): it starts at 8 x 10^ K 
at th e coronal base, peaks a t about 3 x 10^ K at 3 Rq and then falls off with distance as 
r~°'^ (jCasalbuoni et al.lll999l ). The wind speed profile follows from the wind equations with 
given temperature and fiux tube expansion, of the form A{r) = f{r)r'^, with / a function 
whi ch has a maximum clo se to the coronal base and te nds to a finite value at large distances 
(see lKopp fc Holzerlll976l . and lMunro &: Jacksonlll977l ). The same functional form is chosen 



for the expansion in the static part of the atmosphere but different parameters are selected 
in order to obtain realistic values for the magnetic field and its continuos variation at the 
transition region. In the photosphere and chromosphere the profile for Alfven speed is ob- 
tained from the magnetic fiux conservation, B = BqAq/ A{r), and the density profile imposed. 



Following iDmitruk et al.l (l2001bl ) we choose the following model for the nonlinear terms 
in eq. [2] 



NL 



± 



(2) 



^'L(r) 

where L represents an integral turbulent dissipation length and \Z'^\ stands for the total 
amplitude of the Elsasser field integrated over the whole spectrum {Vt) at the point r, hence 

\Z^\ = y//J|zn^)IV^ du]. 

This choice overestimates the transfer rate between high-frequency modes, for which the 
Alfven effect is important (shebalin et al?). In reality the predominant interaction, as will 
be seen below, concerns the lowest frequency reflected mode and the full outward propagating 
spectrum, for which the resonance effects are not important. 

The energy distribution among the modes influences the dissipation rate of all the waves 
coupled. In particular, at a flxed total rms energy, dissipation is reduced if the energy of 
the higher frequency waves is comparable to the lower frequency ones (flatter spectra) with 
respect to the ca se in which rnost of the energy is contained in the low frequency modes 
(steeper spectra) (IVerdini et al.ll2005l ). 

The eqs. [2] can be simplifled including the systematic variation of the Elsasser amplitude 



in a new normalized variable 



N 



{Ma ± l)/v^ 



which re duces to z^r 



the limit of small alfvenic mach number Ma = U jVa —>■ (cfr. iHeinemann fc Olbertlll980l ). 



±z^p!^ in 
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After Fourier transforming in time the hnear equations and adding the phenomenological 
nonhnear term one obtains, 

(U ± K)z^' - tujz% -liU± K)^z^ = -^z± (3) 

(the prime indicates a derivative with respect to r). The numerically integrated equations 
are: 

for the corona, while for the photosphere and the chromosphere one gets: 

The second, third and last coefficient in eqs. HHS] represent the propagation (P), reflection 
(R) and nonlinear dissipation (NL) coefficients respectively (inverse of parallel wavelength, 
reflection scale height, nonlinear length scale). The dissipative feature of the nonlinear 
terms can be shown multiplying the above eq. [3l in its old variables form, by the complex 
conjugate z^* to obtain the evolution equations for the Elsasser energies at a given frequency 
= i|z''=(co')p. On the RHS one gets — |z=''|^|Z^|/L, which is independent of the phase 
difference between the two fields and involves the total amplitude of the fluctuations (the 
same term appears in the equation for a static atmosphere). In the presence of a wind, 
energy flux as conserved quantity is replaced, for linearly propagating waves, by the total 
wave action flux, which may be written as the difference between an outgoing and ingoing 
flux: 

s* = s+ -s- = 



\pUA 



^ + ^"(U + V.)|z-P-^(U-V, 



a/ 



UVa ' ' UVa 

[+, — refer to outward/inward direction and S is the wave action). 



(6) 



The inward wave action density vanishes at the Alfven critical point (Xa ~ 13 Rq, 
where the Alfven speed equals the wind speed), so one may write S* = Sq — Sq = , 
where the index c stands for the critical point, while the index refers to the base of the 
layer. Amphtude and the phase of the outward propagating Elsasser field (z+) at Xa define 
the natural boundary conditions, since the critical point is a regular singular point for the 
incoming wave equation, because phase velocity of the mode vanishes there: total wave 
action density is imposed and the amplitude and phase of z~ can be derived demanding 
the regularity of the solutions at Xa. However, boundary conditions are chosen to assure 
an amplitude of the rms velocity field fluctuations (i.e. summed over the whole spectrum) 
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of ~ 40 km at 1 Rq, as constrained by observations (IBanerjee et al.lll998l ). with an 
assigned spectral distribution: this requires some trial and error since nonlinearity does not 
allow rescaling of the photospheric amplitude by simply rescaling values at the critical point 
Xa- The shape of photospheric spectrum is imposed approximately thanks to the quasi-linear 
properties of the waves in the photosphere-chromosphere layer (small wave amplitudes) and 
the fact that transmission and nonlinearity yield frequency-independent evolution in the 
low corona, as shown in the next section. Given a slope p at the Alfvenic crit ial point, the 
transmission coefficient of the static layer T{lj) (see 
on it), 

St pcVac |Z, 



Krogulec &: Musielakll 19981 for discussion 



■'Nc\ 



PqVo 



(7) 



aO I'^O 



can therefore be used to correct the initial spectrum 

Ml 



|z+(c.)| 
= |z 



desired spectrum at the photosphere imposing 

to describe the spectrum 32 modes are chosen in the range of frequency between 10 
and 10~^ Hz with increasing resolution at higher frequencies. 



z+(tt;o)| X {uj/ujoy to the 
uJo)\^/T{uj) X {io/iOoy. In order 

Hz 



The phenomenological turbulent length scale varies as L(r) = Lq x A(r), where 
Lq = 34, 000 km is imposed at the coronal base and corresponds to the average size of the 
super granules. The waves are propagated from the Alfvenic critical point forward (to the 
Earth orbit) and backward (till the base of the corona) by integration of eqs. HI The con- 
servation of the energy flux across the transition region allows one to determine the Elsasser 
fields below the discontinuity which are propagated back to the base of the photosphere 
using eqs. [51 



3. RESULTS 
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Fig. 2. — Comparison of the reflection (R, thick sohd hne) and nonhnear coefficient (NL, 
thick dashed hne) normahzed to the propagation coefficient (P) for the outgoing and ingoing 
wave at three different frequencies (10~® Hz, 10~^ Hz, 10~^ Hz; black, blue and green lines 
respectively). Also shown in thin solid line is the contribution of each frequency wave to 
the nonlinear coefficient. The photospheric frequency spectrum is fiat and the boundary 
conditions at Xa are set to get a rms velocity field fluctuation at the coronal base 5u = 
40 km s~^ 
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Following IVellil Il993l we compare the characteristic lengthscales of eqs. HHS] in the two 
layers. First consider the thick lines in fig. [2] which represent the refiection and nonlinear 
coefficients (solid and dashed line respectively) normalized to the propagation coefficient 
for uj = 10~^ Hz, 10~^ Hz, 10~^ Hz (black, blue and green lines respectively) for a fiat 
photospheric spectrum. Refiection has a maximum at the transition region and it falls off 
by a factor of about 100 in the corona (because of the density drop). The zeros in the 
refiection coefficient appearing for both the and the z~ depend on the fact that = 
(approximately in the corona), while the one located at Xa appears only for the backward 
propagating waves since the propagation coefficient becomes infinite there (see eq. H]). 
For the outward propagating wave (left panel) refiection is generally much greater (a factor 
100) than dissipation in the photosphere-chromosphere and in the very low corona (below ^ 
1.2Rq). Further out the nonlinear dissipation is smaller than refiection but of the same order 
of magnitude. For the inward propagating wave (right panel), again refiection dominates in 
the photosphere-chromosphere (by a factor of 10), but in the corona the dissipative coefficient 
is comparable or much greater than the refiection coefficient. 

The relative dissipation of the linearly conserved quantities, as defined below in eq. [HI 
has hence different features in the two layers. In fig. [3] we plot the total wave action density 
for the corona (main panel) and the total wave energy fiux for the static layer (sub panel) 
normalized to their base value for all the frequencies which form the spectrum, i.e. 



S*{r,Lu) 




2 _ 




2 






2 _ 




2 



1 1 r / \Zo\ , +p \ZM , -|2 



with the normalization used to derive eqs. [3l the coefficients appearing in the integral are 
the nonlinear frequency integrated coefficients discussed above. 
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Fig. 3. — Normalized wave action density for the corona as function of distance for 5 
frequencies (10^^ Hz solid line, 10^^ Hz dotted line, 10^^ Hz dashed line, 10~^ Hz dotted- 
dashed line, 10^^ Hz triple-dotted-dashed line). The wave energy flux for the photosphere- 
chromosphere is plotted in the subpanel whit the same line coding. 
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Fig. 4. — Heating rate per unit mass integrated over the whole spectrum as function of 
hehocentric distance. The contribution of the ingoing (dashed hne) and outgoing (dotted 
hne) heating rate is also shown. The transition region (T.R.) and the Alfvenic critical point 
(Xa ~ 13 Rq) are indicated on the x axis. 
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In the upper chromosphere the flux tube expansion is very rapid and reflection is strong, 
both the ingoing and outgoing wave contribute to the damping of the energy flux (compara- 
ble -less than one order of magnitude difference- nonlinear coefficient and wave amplitudes) 
and the relative dissipation is very high. Low frequency modes (lower plot in the subfigure) 
are the most damped (the most reflected) while high frequency modes (higher plots) are the 
less damped. In fact, inspection of eq. [8] reveals that the relative dissipation is quadratic 
in the frequency dependent wave amplitudes (|z^p) which in turn increase with decreasing 
frequency because of the different reflection rate. In the corona, instead, beyond 2Rq the 
dissipation coefficient for the outgoing waves is weaker and their amplitudes grow, reflection 
is weaker as well, and an imbalance between outgoing and ingoing fluxes holds. Only the 
former contribute to the wave action density dissipation since now the dominant quadratic 
dependence in eq. E] con ies form the outgoing mode (see fo r example the approximate con- 



servation form used by ICranmer fc van Ballegooijenl |2005| ) . Note that for all frequencies 



the wave action density decreases at approximately the same rate. It turns out that the 
amplitude evolution is driven mainly by the nonlinear, frequency independent, term in the 
corona and by the reflection, frequency dependent, term in the photosphere-chromosphere, 
a feature we will find again studying the power spectrum evolution. 

In comparison the heating rate per unit mass, an absolute measure of energy dissipation, 
integrated over the spectrum, 

Q _ Q+ O- _ \Z+\^\Z-\ + \Z-\^\Z+\ 
— \ rr~\ ' 

P P P L{r) 

is generally higher in the corona than in the photosphere-chromosphere, as shown in fig. |H 
In the latter layer both the ingoing and outgoing wave contributes to the total amount of 
heating rate, while in the former most of the dissipation comes from the outgoing mode. 
The absolute dissipation is quadratic in the frequency integrated wave amphtudes and in 
the corona outgoing wave are allowed to grow almost undamped (low relative dissipation) 
but the existence of a small seed of ingoing wave assures a large absolute dissipation. This 
is not true in the photosphere-chromosphere, before the rapid expansion of the flux tube, 
where the wave amplitude is small and there is a small imbalance between outgoing and 
ingoing propagating wave amplitudes. 




Fig. 5. — Same as fig. [2] for a pliotosplieric Kolmogorov-like frequency spectrum 
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The effect of a different slope of tlie initial spectrum can be understood analyzing the 
contribution of each frequency to the nonlinear coefficient, plotted in thin lines in fig. [2j 
Starting with a fiat photospheric frequency spectrum results in an approximately equal con- 
tribution to the total nonlinear term in the whole atmosphere, except for the outer corona 
where the nonlinear coefficient for the outward propagating wave is made up of essentially 
backward propagating waves at low frequencies. Note also that the frequency decomposed 
nonlinear coefficient is approximately the same for outgoing and ingoing propagating waves 
in the photosphere-chromosphere since reflection is high enough compared to dissipation. 
It follows that if a Kolmogorov-like photospheric spectrum (E / lu oc lu^^^^) is imposed, the 
nonlinear term is mainly made up of low frequency waves for both counterpropagating waves, 
in both the layers. This is can be seen in fig. [5] comparing the dashed thick lines and the 
solid thin lines: for u > 10~^ Hz the contribution to the nonlinear coefficient is generally 
less then 10%. An exception is found below 2Rq for the outgoing mode, since reflection is 
high even for intermediate frequency wave (see fig. [7| for the photospheric layer) . Note that a 
deap in the (frequency integrated) nonlinear coefficient for the outgoing mode appears below 
the location of vanishing reflection in both the photosphere and low corona, since the energy 
resides mainly in the low frequency mode. 
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Fig. 6. — Compensated power spectrum as function of heliocentric distance for a photospheric 
kolmogorv-iike initial spectrum. Each curve is labelled with the corresponding heliocentric 
distance in unit of Rq. Prom top to bottom, solid lines indicate the photospheric base, the 
T.R., Xa and 1 AU. 
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This separate behavior in the two layers has strong consequences on spectral evolution. 
In fig. [6] the (compensated) total power in the fluctuations is plotted for different heliocentric 
distances. An almost Kolmogorow-like spectrum is imposed at the base of the photosphere 
with a procedure described at the end of section [21 At very low frequencies the spectrum 
practically does not evolve, in the whole domain, while there is a tendency to steepen at 
low-intermediate frequencies (10~^ Hz < a; < 10~^ Hz). The behavior at high frequencies 
is quite complicated. Some irregularities appear very close to the base of the photosphere 
and the overall tendency is that of flattening. Note, however, that most of the changes 
in the shape occur in the photosphere-chromosphere where the waves display a strong fre- 
quency dependent behavior. This makes the spectral evolution very similar to the linear 
case, below the transition region (except the energy level of the spectrum), and the appear- 
ance of the irregularities can be interpreted by means of the linear analysis. Accordingly, 
in flg. [7] we plot the transmission coefficient, defined in eq. [71 as function of frequency for 
the photosphere-chromosphere. Note that the transmission is constant at low frequency, 
decreases at intermediate frequencies and increase again at high frequencies, where several 
transmission peaks appear: basically all spectral evolution is qualitatively reproduced. 
The pe aks origina te form the discontinuity in the refiection scale height at the transition 
region (1Vellilll993l ). In fact the amplitude of the refiected waves shows some nodes inside the 
domain and when their location coincides with the base of the photosphere t he transmissio n 
is enhanced (a condition which depends on the frequency of the waves, see iHoUwegl Il978l ) . 
When nonlinearities are in troduced the locatio n of the nodes depends also on the wave am- 



plitude imposed at Xa (see lVerdini et al.ll2005l ) and similarly if these nodes are located near 



the base of the photosphere the irregularities in the spectrum appear. 



The slope of the spectrum imposed at the photosphere has negligible effects on the total 
power spectral evolution, however it changes the amount of energy residing in the ingoing 
and outgoing mode (or in the kinetic or magnetic fluctuations) at large distances and some 
constraints on the slope can be obtained using the available observational data. In flg. [8] the 
Elsasser energies E"^ integrated over the frequency spectrun a are plotted (solid and dashed 
line respectively) along with the Ulysses and Helios data (iBavassano et al.l l2000bl ) . for a 
Kolmogorov (thick lines) and a flat (thin lines) initial slope with 6u = 40 km s~^ at the 
coronal base. Both the data and the expected slopes are reproduced by the Kolmogorov-like 
photospheric spectrum while the flat one has to high outgoing energy and too low ingoing 
energy. The effect of high energy at high frequency waves is that of dissipating the inward 
wave, since they are mainly outward propagating: as a result outgoing waves are allowed 
to propagate almost undamped and their energy content is hence higher. Note that in the 
Kolmogorov case a deap, very close to the coronal base, appears as a signature of vanishing 
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ingoing waves, a feature of the low frequency reflected waves. This results in a vanishing 
absolute dissipation (heating) which is not found for the flat case and has important conse- 
quences for the acceleration and heating of the solar wind. 




Fig. 8. — Frequency integrated Elsasser energies as function of heliocentric distance for a 
photospheric Kolmogorov spectrum with 5u — 40 km s~^ at the coronal base. Symbols 
indicate observational constraints (see text for explanation). 
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In the following we consider only a Kolmogorov spectrum. In fig. [9] the root mean 
square amplitude of velocity field fiuctuation integrated over the whole spectrum is plot- 
ted as function of heliocentric distance (solid line, in dotted line we plot also the mag- 
netic field fiuctuation in velocity unit) along with some observational data (taken from 



Cranmer fc van Ballegooiienll2005l . to which we address for comment on the data set): 



Filled di amonds are nontherm al line broadening velocities measured by SUMER on 
the disk Jwilhelm et al.lll995h . 



Crosses are nonthermal velocities derived from SUMER observations above the solar 



limb (IBanerjee et al.lll998l ) 



The box represents the upper and lower limit given by lEsser et al.l Il999l from UVCS 
off-limb data 



Stars are early measurements from [Armstrong fc Wodll981 



The bars are recent measureme nts of transverse velocity field fiuctuation using radio 
scintillation (jCanals et al.ll2002l ) 



Filled bars are the Helios and Ulysses data for the Elsasser energies, from lBavassano et al. 



2000bl . rewritten in term of the velocity field fiuctuation assuming equipartition between 



magnetic and kinetic energy. 



Note that the Helios and Ulysses data are obtained averaging over periods >1 hour (corre- 
spondig to < 10~^ Hz) while all the other points in the figure refers to rms values. 




Fig. 9. — Root mean square amplitude Su and Sb (in velocity units) as function of heliocentric 
distance for a photospheric Kolmogorov spectrum with 5u — 40 km s~^ at the coronal base. 
Symbols indicates observational constraints (see text for explanation). 
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The overall agreement is quite good, even if data suggest a smaller power (more dis- 
sipation) just above the T.R. and more power (less dissipation) at about 2Rq. Note that 
because of the equipartition assumption the Helios and Ulysses data disagree with the in- 
tegrated quantities (the correct comparison has already been made above in fig. [H]). As 
noted by lCranmer fc van Ballegooiienll2005l the longitudinal velocity fluctuation data (filled 
diamonds) agree very well with the magnetic field fluctuation amplitudes (dashed line) and 
could indicate wave coupling among transverse and longitudinal mode. At leading order 
compressive effects are driven by the magnetic pressure originating from the incompressible 
fluctuations and represent a way for Alfven waves to get rid of the energy excess above the 
T.R.. If these compressional waves are isotropic and suffer some dumping via shock forma- 
tion, or other processes active in the low corona, they can reproduced the measured parallel 
5u and supply the heating needed by current model of wind acceleration. 



4. CONCLUSIONS 

In this paper we have modeled the nonlinear evolution of Alfven waves propagating 
through the photosphere, the corona and the solar wind till 1 AU. Nonlinear interactions 
occur between outward propagating and reflected waves, and it is assumed that a nonlinear 
cascade develops preferentially in a direction perpendicular to that of propagation, which we 
take to coincide with the direction of the mean radial magnetic field. 

While the phenomenological nonlinear term acts as a dissipative sink for both outward and 
inward waves, independently of the wave frequency, reflection, provided by the stratification 
of the layer, is generally strong at low frequencies and decreases with increasing frequency. 



We find that most of the heating occurs in the low corona (below the Alfvenic criti- 
cal point), while very little power is dissipated below the transition region. For reasonable 
velocity field fluctuations at the base of the photosphere a sufficient amount of energy flux 
is transmitted through the transition region. The adopted frequency coupling is not able 
to reproduce the observed spectral slope and evolution in the Alfvenic range even though 
frequency integrated data at large distances constrain the outer spectrum to be steep (-5/3 
slope). The modification of the frequency spectrum occurs mainly in the chromosphere and 
in the photosphere since waves experience a strong reflection at all the frequencies consid- 
ered, while in the corona and the solar wind the spectrum maintains approximately the same 
shape one finds at the coronal base. 
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Nonlinear dissipation based on reflection acts in different ways depending botfi on tfie 
(ingoing and outgoing) wave amplitude and on the layer considered. In the corona, reflec- 
tion is not very high but the outgoing wave amplitude is allowed to grow, so that the wave 
evolution is driven by the nonlinear interactions (all the modes evolve in the same way) and 
one finds a strong heating rate in the Sub-Alfvenic corona. 

In the photosphere-chromosphere a strong reflection rate, combined with small wave ampli- 
tudes, leads to an evolution similar to the linear case, which depends on frequency, and a 
small heating rate. 

As a result most of the wave energy dissipation takes place in the first 4 solar radii above the 
coronal base. The driving modes for dissipation are the modes which experience the biggest 
refiection, generally low frequency modes. However depending on the model of atmosphere, 
i.e. on its characteristic scale height, and on the energy distribution, i.e. fiat or steep spectra, 
intermediate frequency modes can be important as well. 

The spectral shape varies mainly below the transition region, it steepens at low-intermediate 
frequencies (10^^ Hz < < 10^'^ Hz) and fiattens at high frequencies {uj > 10^^ Hz) show- 
ing the characteristic features (energy peaks and frequency distribution) one finds in the 
transmission coefficient (linear behavior). In the corona instead, it maintains approximately 
the shape one finds beyond the T.R., because of the form of the nonfinear term adopted. 
The very low frequency range (a; < 10~^ Hz) practically does not evolve in the whole layer 
and it keeps the original slope at the photosphere. With this model of nonfinearities one can 
conclude that the spectrum one finds at 1 AU is basically the same spectrum at the base of 
the corona. 

The input spectrum at the photosphere, whatever the shape is, is instead strongly modified 
by the transmission properties of the atmosphere below the transition region (independently 
of the model used for the nonlinear interaction). The energy peaks in the spectrum, result- 
ing from an enhanced transmission at high frequencies, indicate that, even in presence of 
nonlinar interactions, the photospheric layer act as a filter for the energy injected through 
photospheric footpoint motion, if a smoothing of the forcing frequency is to be present, it 
must occur in this highly stratified layer. 

The data at large distances suggest that energy at high frequency should be very low, however 
we find an energy increase at high frequency. Since the spectral evolution in corona depends 
also on the approximate frequency coupling contained in the nonlinear term, constraints on 
the photospheric input spectrum can not be given safely. Given that high frequency waves 
are transmitted through the T.R. and are quite energetic in the very low corona some other 
mechanism must be invoked to dissipate high frequency waves or a better modeling of the 
nonfinearities, which we plan to do in future works. 
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As first pointed out by iHollweallQSll . such high frequency energy reservoir can be the source 
for plasma heating processes operating in the low corona. Note that not only the peaks 
contribute to the energy budget, but the general flattening of the spectrum is important as 
well. 

A comparison with measurements of 5u suggests that the model can be considered a very 
good approximation in the outer corona and solar wind while, despite the good agreement 
found in the low corona, some other processes must be invoked to reproduced the observed 
features below the alfvenic critical point, such as compressible effects and wave coupling, 
especially in the chromosphere and photosphere. Other models of turbulent transport have 
been constructed t o fit the decay of turbulence with d istance from the sun in the solar 
wind beyond 1 AU fISmith et al.ll200ll : iBreech et al.ll2005l ) as well as to explain the extended 
heating in this region. Here the Alfven speed can be neglected in the transport of the fluc- 
tuations, so that in some sense our model equations should be consistent with theirs, when 
rewritten in terms of the second order moments. A generalisation of turbulence transport 
equations, consistent both in the corona, acceleration region and solar wind is a topic of 
current research. 
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